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Abstract 

We explain how an analytic continuation we have developed recently is applied to determine 
the residues of the nucleon resonance poles within a dynamical coupled-channel model of meson- 
baron reactions. A procedure for evaluating the electromagnetic N-N* transition form factors at 
resonance poles is developed. Illustrative results of the obtained N* — ► ttN, jN transition form 
factors for Pn, P33, and D13 nucleon resonances are presented and compared with previous results. 
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I. INTRODUCTION 



In a recent paper [lj, we have developed an analytic continuation method to determine 
nucleon resonances within a dynamical coupled- channel model of meson-baryon reactions 
(MSL). The method has been applied^ to extract 14 nucleon resonances from the ttN model 
developed in Ref.[H (JLMS) which has been extended to investigate ttN^N — > 7T7r./v{5j], 
7 iV -> ttNQ and N(e,e'ir)NQ reactions. The purpose of this paper is to explain how the 
residues of the extracted resonance poles are determined from the predicted ttN — > ttN and 
7*iV — > ttN amplitudes. 

In section II, we briefly review the analytic continuation method developed in Ref.[l[. 
Section III is devoted to explaining the determination of the residues of the nucleon resonance 
poles. Illustrative results for Pn, P 33 , and _Di 3 nucleon resonances are presented in section 
IV and compared with the results from other analysis. A summary is given in section V. 



II. ANALYTIC CONTINUATION METHOD 



Within the MSL formulation[2j, the partial wave amplitudes of two-body meson-baryon 
reactions can be written as 

7>, a (p',p;£0 = tj3 >a (p',p;E)+tX a (p',p;E), (1) 

where a, (3 represent the meson-baryon (MB) states 7 AT, ttN, rjN, pN, aN, 7rA, and 

t* a (p',p;E) = ^pAp';E)[G N ,(E)} hl T a ,(p;E) (2) 

with 

[G N \] hl (E) — (E — m N: )5 hJ - E i}j (E) . (3) 

Here i, j denote the bare A^* states defined in the Hamiltonian. rriN* are their masses. 
The first term (called meson-exchange amplitude from nowon) in Eq.(TjQ) are defined by the 
following equation 

t/3 >a (p',p;E) = vp !a (p',p)+ / dqq 2 J2 v P,7(p' ',Q', E )Gj(q,E)t J>a (q,p;E) (4) 

where vp tCt is defined by the meson-exchange mechanisms, and G 1 (q,E) is the propagator 
for channel 7. The dressed vertices and the energy shifts of the second term in Eqs.(T5])-([3]) 
are defined by 

T aJ (p;E) = T aJ (p)+ / dqq 2 J2 t an(p\^ E )G^(q,E)T^ J (q) (5) 
£(£) M = / dqq^^MG^E^q). (6) 

J C ry 



where T aji (p) defines the coupling of the z-th bare A^* state to channel a. 

To search for nucleon resonances, we have developed [lj] an analytic continuation method 
to find poles of the scattering amplitude T a ^{p' , p\ E) on the unphysical sheets of the complex 
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FIG. 1: The shift of the on-shell momentum (open circle) of the two-particle Green function Eq. ([7]) 
as energy E moves from a real value above the threshold energy to a complex value with negative 
imaginary part. C[ in (a) or C\ in (b) is the integration path for calculating Eqs.dU-Q amplitude 
for E on the unphysical plane. 



energy E'-plane. For multi-channel reactions considered here, there can be many poles 
associated with a single resonance on different unphysical sheets. The pole nearest to the 
physical sheet is supposed to play a dominant effect on physical observables, and is called 
the resonance pole traditionally. The other poles are called shadow poles. 

Since v a ^ and the bare vertex T a ^ are energy independent within the MSL formulation, 
the analytic structure of the scattering amplitude defined above as a function of E is mainly 
determined by the the Green functions G~,(q, E). Thus the key for selecting the amplitude on 
physical sheet or unphysical sheet is to take an appropriate path of momentum integration 
C in Eqs. ([T]) -([51) according to the locations of the singularities of the meson-baryon Green 
functions G Q (p,E) as E move to complex plane. This can be done independently for each 
meson-baryon channel. For channel with stable particles such as ttN and rjN, the meson- 
baryon Green function is 

Gu B (E, P )= E _ E J ) _ EB{py (7) 
which has a pole at the on-shell momentum po defined by 



E = V m M+Po + V m B+Po- (8) 

As an example, let us consider the analytic continuation of the amplitude to the unphysical 
sheet of the MB channel when the energy E is above the threshold Re(E) > mg + rriM and 
Im(E) < 0. The on-shell momentum po for such a E is on the second and the fourth quadrant 
of the complex momentum plane. As Im(E) becomes more negative as illustrated in Fig. 
[TJ the on-shell momentum (open circle) moves into the fourth quadrant. The amplitude on 
the unphysical sheet can be obtained by deforming the path C into C[ or equivalently C\ so 
that the on-shell momentum does not cross the integration contour. We note here that for 
the energy below threshold (Re(E) < tub + tum) the path G\ will give amplitudes on the 
physical sheet. 
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FIG. 2: Contour C2 for for calculating Eqs.([l|)-([6]) for E on the unphysical plane with the unstable 
particle propagators, such as Eq.([9j) for irA channel. See the text for the explanations of the dashed 
line and the singularity p x . 



For the channels with unstable particle such as the 7rA as an example, the Green function 
is of the following form 

G ^ E ^ = E-EM-E A (p)-z A {E, P y (9) 



where 



c 3 E- EM - [(E„(q) + E N (q))*+p^/i- 

(10) 

The 71 A Green function Eq. (Q has a singularity at momentum p = p x , which satisfies 

E - E„(p x ) - E A {p x ) - Z A (p x , E) = 0. (11) 

Physically, this singularity corresponds to the ttA two-body 'scattering state'. There is also 
discontinuity of the ttA Green function associated with the ttttN cut in Sa, as shown in the 
dashed line in Fig. [21 where po is denned by 

E = K(po) + [Kr + m N ) 2 + pi} 1 ' 2 . (12) 

Therefore, for Re(E) > m B + m M , Im^ + m N , the integration contour C must be chosen to 
be below the titiN cut (dashed line) and the singularity p x , such as the contour C2 shown 
in Fig l2j for calculating amplitudes on the unphysical sheet. 

The singularity go of the integrand of Eq. ffTUj) depends on the spectator momentum p 

E - E n {p) = [(E^qo) + E N (q )) 2 + p 2 ] 1/2 . (13) 

Thus go moves along the dashed curve, illustrated in FigEl when the momentum p varies 
along the path C2 of FigEl To analytically continue Ea(p, E) to the unphysical sheet, the 
contour C3 of Eq. ffTUj) must be below go- A possible contour C3 is the solid curve in Figj3l 
We emphasize here that we can deform the contour C only in the region where the poten- 
tial v a ^{p' \p) and the bare N* vertex T MBtN *(p) are analytic. The contours described above 
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FIG. 3: Contour C3 for calculating the A self energy Eq. lfTU]) on the unphysical sheet. Dashed 
curve is the singularity qo of the propagator in Eq. (|13p , which depends on the spectator momentum 
p on the contour C2 of Figj2j 



only from considering the singularities of MB and 7T7riV Green functions. Thus they must 
be modified according to the analytic structure of the considered v a ^(p',p) and T MB)N *(p). 
Within the MSL formulation, the t-channel meson exchange potential ^Vb' mb(p iP) nas 
singularities at 

A 2 -{P-Pf = (14) 

with A = E M r(p') — E M {p) or Es'ip') — E B (p). The form of r mb,n* (p) is chosen such that 
its singularity is at the pure imaginary momentum. Thus the contours have to be chosen to 
also avoid these singularities. As an example we show in FigJH the singularities associated 
with the 7rA channel at E = 1357 — 76iMeV. The dotted line for irnN cut and the circle 
shows px are the singularities from the Green's function, as discussed above. The most 
relevant singularity of the meson-exchange potential in our investigation of electromagnetic 
pion production amplitude is due to the t-channel pion exchange of 7iV — > 7rA, which is 
shown as the dashed-dot curve. Thus the integration contour has to be modified to the solid 
curve in FigJH 

III. EXTRACTION OF RESONANCE PARAMETERS 

The resonance energy (M# = M — iT/2) is the position of the pole of the scattering 
amplitude which is on the sheet nearest to the physical sheet. In principle the resonance 
pole can be found in the meson-exchange amplitude t and/or resonance amplitude t R of Eq. 
(CQ). Within the nN model developed in Ref. [4j (JLMS), we findjsj that resonance poles 
are only from t R . We therefore will only explain how the residues of resonance poles are 
extracted from this term. 

The poles of t R are found from the zeros of the determinant of A^* propagator defined by 
Eq.© 

A{E) = det[G N \(E)} = 0. (15) 
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FIG. 4: The contour (solid curve) for calculating electromagnetic matrix element, po and p x are 
the singularities shown in Figj2j The dashed-dot curve is the singularity of the pion-exchange 
jN — > irN matrix element at E = 1357 — 7Qi MeV. 



Near the resonance energy Mr, N* Green function can be expressed as 



XiXj 
E-M, 



R 



(16) 



where i,j denote the bare N* state in the Hamiltonian and Xi represents z-th component of 
the dressed N* and satisfies 



J2(G N *(M R )-% Xj = YH(MR-mN*)Si j -X(M R ) ij } Xj = 0. 



(17) 



3 3 

If there is only one bare N* state, it is easy to see that 

1 

X = 



where E'(Mr) = [d'E/dE]E=M R - If we have two bare N* states, we find that 



(18) 



Xi 



X2 



M R - m N ; - T, 22 {M R ) 



Si2(M fl ) 



Mr - m N * , -E 22 (Mr 



■Xi 



(19) 
(20) 



where A' (Mr) = [dA/dE] E=MR can be evaluated using Eq. ffl~5l) . 

We now examine how the residues Xi can be used to see the effects of resonance poles on 
the full amplitude T a ^ defined by Eq.([T]). Near the resonance pole we can perform Laurent 
expansion of the on-shell amplitude 



Ta,p(Pa,P% M R 



-pRrR 

1 a L f3 



E-M 



R 



+ c°+cUe-m r ) + .... 



(21) 
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For the cases that the resonance poles are from t R term of Eq.flT]), one can see from the 
definitions Eqs.©-© and Eq. (jT6l) that the vertex functions in the residue of Eq. (l2T]) is 
determined by the dressed vertex T a j defined by Eq.fJS]) 

ff = ExAiba,^)- (22) 
3 

The terms C^p and C\ a in Eq. f[2~Tj) also depend on the matrix elements of meson-exchange 
amplitude t of Eq.Q, but will not be discussed here. 

The residue of the nN elastic scattering amplitude F^n^n characterizes the strength of 
the coupling of the resonance with nN channel. Using the standard notation, we have at 
near resonance position Mr 

F nN ,ME) - s " [ M^E 1e ^ Mr ' (23) 

where S^n^n is the partial-wave S-matrix. In terms of the normalization of JLMS model, 
we have 

F-kN^nKMr) = -7r — T wNi7rN (po,p ,M R ) (24) 

Mr 



We thus have 



fti# = ,*(^w teV (25) 



The elasticity of a resonance is then defined as 

* = (26) 

The electromagnetic N-N* transition form factor is defined by matrix element of the 
electromagnetic currents between nucleon and iV* 

A 3/2 (Q 2 ) = X <N*,s z = 3/2\- J(Q 2 ) -e +1 \N,s N = 1/2 >, (27) 

~>2\ v ^ t\t* „ 1 /oi 7/ rft\ 



A 1/2 (Q 2 ) = X < N*,s z = l/2\- J(Q 2 )-e +1 \N,s N = -1/2 > ,, (28) 

1/2 

where e+i = — (x + iy)/V2, 



S 1/2 (Q 2 ) = X < N*,s z = l/2\J°(Q 2 )\N,s N = 1/2 > , (29) 



X 



l E N (q) 1 



V % V2K 

with K = (M R — m 2 N )/{2MR). The above definition was originally introduced for the 
constituent quark model [8]. If < iV*| is a resonance state, then the above expression at 
resonance position Mr must be evaluated by using Eq. f[2~2"j) . We thus have 

A 3/2 (Q 2 ) = XX'^Xj^NjiQ^M^X^ = 1,X N = -1/2) , (30) 

(31) 

where 



X' 



, (2j + l)(27r)3(2g ) 



V 47T 

Here the additional factor X' is due to our normalization of the vertex function T. The above 
helicity amplitudes are in general complex number. Au2 and S1/2 have similar expressions. 
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IV. ILLUSTRATIVE RESULTS AND DISCUSSIONS 



In this section, we illustrate our procedures by presenting the results for the pronounced 
resonances in P 33 , D± 3 and in the most complex P n partial waves. Their pole positions were 
determined in Ref. [3| and are listed in Table [H Our results for P 33 (1232) and Pi 3 (1520) 
agree well with the values listed by Particle Data Group [9] (PDG). For P n channel we 
found three poles below 2 GeV. Two of them near 1360 MeV are close to the tc A threshold. 



This finding is consistent with the earlier analysis of VPI [10(| and Cutkosky and Wang [11 



and the recent analysis by the GWU/VPI[12| and Juelich[13j groups. We have shown in 
Ref. [jl] that all of the three Pn resonances listed in Table [TJ correspond to a single bare 
iV* state at 1736 MeV. This is a dynamical verification of the resonance pole-shadow pole 
relation in coupled-channels reactions, as discussed by Eden and Taylor [14j, Kato[l5|, and 



Morgan and Pennington [16 



The extracted residues Re 1 ^, defined in Eq. ([2"51) . for ttN amplitude are compared with 
some of the previous works in Table [UJ We see that the agreement in P 33 and Pi 3 are 
excellent. For Pn, there are significant differences between four analysis. As discussed in 



Ref. [11[ , it could be mainly due to the differences in the employed reaction models. On the 
other hand, the difference between the predicted Pn amplitudes at W > about 1.6 GeV 
could be the reason why the third Pn pole is not found in Juelich analysis. 

From the values of Mr of Table [TJ and R of Table [Til we can evaluate the elasticities r) e 
using Eq. fTSB"]) . The results are also listed in Table [TJ We see that our results agree well with 
PDG values, while some investigations are needed to understand better the comparisons for 
the two Pn poles near 1360 MeV which are close to 7rA threshold. 

To extract helicity amplitudes using Eq. (j3U|). we use the multipole amplitudes calculated 
from using the parameters determined in Ref.[7|. Our results at photon point are listed in 
Table ITTT1 We observe that the re al p arts of our results for P 33 and Di 3 are in good agreement 
with several previous results 17], [l8|, [l9|, |20j. The large differences in Pn indicate that more 
investigations are needed to understand the differences between our resonance extraction 
method within a coupled-channel model and other methods which are mainly based on the 
Briet-Wigner parametrization of single channel K-matrix amplitudes. 

For P 33 we can use the standard relation [2l[ to evaluate the N-A magnetic transition 
form factor G* M in terms of helicity amplitudes. The real parts of our results are the solid 
circles in FigJ^l which are in good agreement with the previous analysis. In the same figure, 
we also show that the imaginary parts of our results are much weaker. This result and 
the results of Table II II I suggest that we can only make meaningful comparisons with the 
results from analysis based on the Briet-Wigner parametrization of single channel K-matrix 
amplitudes only for the cases that the imaginary parts are small. This turns out to be also 
the case of the Pi 3 (1521) resonance. In Fig[6], we see that the real parts of our A 3 / 2 and Au 2 



are in good agreement with the results from CLAS collaboration [23J. The large differences 



in 5*1/2 perhaps are mainly from the fact that the longitudinal parts of the amplitudes can 
not be well determined with the available data. 

For P 11; the imaginary parts of the calculated helicity amplitudes for the three poles 
listed in Table [TJ are very large. Thus it is not clear how to compare our results with 
previous results. We thus show both the real parts (solid circles) and imaginary parts (solid 
triangles) in Fig. [71 It seems that the structure of A*(1356) and A*(1364) are similar. In 
particular their real parts of A1/2 chan ge s ign at low Q 2 , similar to what haven been seen 
in the results from CLAS collaboration [23]. However, because of the double pole structure 
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TABLE I: Resonance poles (ReM R , -ImM R ) MeV and elasticity rj e (Eq.([26])) extracted in Ref.@] 





M R (EBAC-DCC) 


Mr (PDG) 


r/ e (EBAC-DCC) 


% (PDG) 


P33 


(1211, 50) 


(1209 - 1211, 49 - 51) 


100% 


100 % 


D 13 


(1521, 58) 


(1505 - 1515 , 52 - 60) 


65 % 


55 - 65 % 


Pn 


(1357, 76) 
(1364,105) 
(1820, 248), 


(1350 - 1380, 80 - 110) 
(1670 - 1770, 40 - 190) 


49 % 
61 % 

8 % 


60 - 70 % 
10 - 20 % 



TABLE II: The extracted ttN residues Re irt> denned by Eq.(l25l 





EBAC-DCC 




GWU-VPIfl2] 




Cutkosky[ll] 




Juelich[13] 






R 


<\> 


R 


4> 


R 


<t> 


R 


4> 


P 33 (1210) 


52 


-46 


52 


-47 


53 


-47 


47 


-37 


D 13 (1521) 


38 


7 


38 


-6 


35 


-12 


32 


-18 


Pn(1356) 


37 


-111 


38 


-98 


52 


-100 


48 


-64 


(1364) 


64 


-99 


86 


-46 










(1820) 


20 


-168 






9 


-167 







and the large imaginary parts, more detailed investigations are need to make meaningful 
comparison with previous results. 



V. SUMMARY 

In this paper, we have briefly reviewed the analytic continuation method developed in 
Ref. [lj and explained how it is used to determine the residues of nucleon resonance poles. 
To illustrate our method, we have presented the results for resonances in P 33 , P>i 3 , and P\\ 
partial waves. 

For residues associated with irN channel, we agree with most of the previous results 0, 
f o r P 3 3(1232), £>i 3 (1521) and two P n poles near 1360 MeV. For Pn(1820), the 



calculated elasticity ~ 8% agree well with the value ~ 10 — 20% of PDG and Ref. [ill], while 
this resonance is not reported in the analysis of Refs. fiol. 12. 13l . 



For residues associated with 7^ channel, the corresponding helicity amplitudes for P 33 
(1232) and _Di 3 (1521) are dominated by their real parts which are in good agreement with 
other analysis based on the Briet-Wigner parametrization of K-matrix amplitudes. For 
Pn resonances, the extracted helicities amplitudes have large imaginary parts and more 
investigations are needed to compare our results with previous analysis. 

Our next necessary task is to examine how to define the residues associated with unstable 
7rA, pN, and aN channels. Our effort in this direction along with our complete results for 
the 14 nucleon resonances extracted in Ref.0] will be reported elsewhere. 
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